Coexistence of superconductivity and incommensurate magnetic order 
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The influence of incommensurate spin density waves (SDW) on superconductivity in unconven- 
tional superconductors is studied by means of the Bogolubov-de Gennes (BdG) equations. Ex- 
ploiting translational symmetries of a magnetically ordered two-dimensional system we propose an 
approach that allows to solve the BdG equations on much larger clusters than it is usually possible 
for inhomogeneous systems. Applying this approach we demonstrate that the presence of incommen- 
surate spin density waves induces real-space inhomogeneity of the superconducting order parameter 
even in the absence of external magnetic field. In this case a homogeneous order parameter of 
the Bardeen-Cooper-Schrieffer-type superconducting state is slightly modulated, or equivalent ly, a 
small fraction of the charge carriers form Cooper pairs with non-zero total momentum. However, 
when a suflficiently strong magnetic field is applied, the homogeneous component of the order param- 
eter is suppressed and the system transits to the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, 
where the order parameter oscillates changing sign. We show that for s-wave pairing the presence of 
external magnetic field diminishes the destructive infiuence of the SDW order on superconductivity. 
A simple explanation of this effect is also proposed. 

PACS numbers: 74.25. Ha, 74.70. Tx 



I. INTRODUCTION 

The interplay of magnetism and superconductivity has 
been investigated both theoretically and experimentally 
for many decades. In 1957 Ginzburg demonstrated that 
superconductivity and long range ferromagnetic order 
compete with each other making their coexistence al- 
most impossible^ However, it turned out that this in- 
compatibility can be overcome by spatial separation as 
it takes place, e.g., in some Chevrel phases. The situa- 
tion is particularly interesting in the case of unconven- 
tional superconductors, where the superconducting phase 
is often close to a magnetic one and, what is even more 
important, magnetic and superconducting phases are 
formed by the same electrons. Among others, unconven- 
tional superconductors include cuprates, heavy fermion 
systems^^ and the recently discovered iron pnictides^^ — 
A microscopic description of the interplay between mag- 
netism and superconductivity still represents a complex 
and unsolved problem. 

There are, however, unconventional superconductors 
with coexisting magnetic and superconducting orders. 
Such a coexistence is observed, e.g., in iron pnictides^ and 
in heavy fermion CeRhlns.^ In these systems crossovers 
from a purely magnetic phase to a phase with coexisting 
orders and then to a purely superconducting phase (or 
in the opposite direction, depending on the control pa- 
rameter) are observed. These examples show that mag- 
netism can coexist with superconductivity, but they do 
not say much about the mutual interplay of these or- 
ders, i.e., whether they coexist because of a cooperation 
or despite a competition. The competition can be de- 
duced indirectly from the phase diagrams (see e.g., Ref. 
0) where it shows up as a negative correlation between 
the magnetic and superconducting transition tempera- 
tures. The onset of magnetic order in the vicinity of 
superconducting vortex cores^— gives additional sup- 



port for this competition. Some additional light can 
be shed on this problem by studying a system sim- 
ilar to CeRhlns but with cobalt instead of rhodium. 
CeCoIns is unusual in that the magnetic order occurs 
only inside the boundaries of the superconducting phase 
and vanishes together with superconductivity at the up- 
per critical field.^ This indicates that in CeCoIns there 
is a mutual cooperation between superconducting and 
magnetic orders. The physical origin of the high-field 
and low-temperature (HFLT) phase with coexisting or- 
ders remains controversial, especially with respect to 
its relation to the FFLO super conduct ivityJ^^ri^ While 
the early results on CeCoIus were interpreted mostly 
in terms of the FFLO scenario^— the discovery of 
the incommensurate SDW in the HFLT regime^i^i^^ — 
raised the question whether this phase is of magnetic or 
FFLO nature. Contradictory conclusions can be found 
even in the very recent theoretical^^ — and experimental 
paper s^^i^ii^ii^i^— The observed field independence of 
the wave-vector associated with SDW has been consid- 
ered as a key argument against the FFLO scenario*^ 

The FFLO state is characterized by the formation of 
the Cooper pairs with nonzero total momentum q and 
the spatially inhomogeneous order parameter (OP) of 
the Fulde-FerrelU^ ^(^) ^ exp{ir • q) or the Larkin- 
Ovchinnikovi^ A{r) ~ cos(r • q) types. Usually, there 
are more than two equivalent gr-vectors which give the 
same upper critical field for the FFLO state, while the 
lowest free energy is obtained for OP being a linear com- 
bination of several plain waves^— with a complicated 
spatial modulation of the OP. Since any additional mod- 
ulation of the OP in the real space requires an additional 
components in the momentum space, inclusion of several 
momenta q should be necessary in the presence of mech- 
anisms which break/modify the translational invariance, 
like impuritie a^^i^^i^^ or incommensurate SDW. 

The spatial structure of the FFLO OP has been de- 
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termined mostly from the solution of the BdG equations 
in the real spacei^^~— Such an approach is best suited 
for investigations of the FFLO phase in the presence of 
impurities or vortices when the translational invariance 
is broken. However, for a coexisting SDW and FFLO 
phases calculations in the momentum space allow one to 
benefit from a possible translational invariance in the di- 
rection perpendicular to the SDW modulation. In this 
paper we develop such an approach and study the in- 
fluence of SDW on s and d-wave superconductivity. By 
carrying out numerical calculations on systems up to 10^ 
sites we show that incommensurate SDW itself favors 
pairing with nonzero momentum for both the symme- 
tries. What is interesting, in the presence of incommen- 
surate SDW, the 5-wave superconductivity persists up 
to much higher magnetic fields than in systems without 
SDW. We present simple arguments explaining this ef- 
fect. 



II. MODEL AND APPROACH 

While the spatial structure of the magnetic order 
in HFLT phase has been determined from several 
experiments ^i^i^i^ the spatial structure of the super- 
conducting OP remains unknown. Therefore, the experi- 
mental data on SDW will be taken as a phenomenological 
input in our calculations. Solving the BdG equations for 
the superconducting OP we determine how the assumed 
SDW affects formation of Cooper pairs with various total 
momenta. We investigate the following Hamiltonian on 
a two-dimensional square lattice 



(1) 



where Hs (Hd) represents pairing interaction responsible 
for 5-wave (d-wave) superconductivity and 



Ho = -tJ2 



cLCi 



^ {<T [/i + M{Ri 



(2) 

Here, cj^ creates an electron with spin a at site i, t is 
the hopping integral between the nearest-neighbor sites 
and /i is the chemical potential. We focus on the role of 
the external magnetic field h and incommensurate SDW 
order MjRj) = Mocos{Ri • Qsdw)- Following Refs. 
la, H, S and im we take either Qsdw = (Q,Q) or 
QsDW = (Q,7r). Hamiltonian (j2j) includes the Zeeman 
pair breaking but it neglects the orbital effects of mag- 
netic field, which are the most effective pair-breaking 
mechanism in many superconductors. Here, however, we 
focus on heavy fermion systems where this mechanism is 
ineffective due to the large electron effective mass. This 
mechanism does not play a role also in layered systems 
provided the field is applied parallel to the planes. 

The mean-field form of the on-site pairing interaction 
for s-wave superconductivity reads 



H. = 



B.C. 



(3) 



with = Vs{ciiCi^) and Vg < 0. In the case of inter-site 
pairing we assume 



where Af = ^(Q+c^^Qt " Ci^a^cn) for a G {x,y} and 
Vd < 0. This form of the interaction Hamiltonian al- 
lows for an arbitrary inter-site singlet pairing (e.g., the 
extended s-wave). However, the experimental results 
suggest d-wave superconductivity and we restrict further 
analysis to this type of pairing. Representing the super- 
conducting order parameters by their Fourier transforms 

A|") = ^A(")exp(iq-Ri), (5) 
q 

one obtains the Hamiltonian in the momentum space 

- ~Y' ('^fc<^'^*'-QsDW,<T + cL'^fc+QsDW,<T) , (6) 



H. = Y.A, cld^^^, + H.c. - - ^ I A,^ (7) 

kq q 

= E E ^"C^' 9) ^U^-k+,i + H-c. 

kq a 

-^Ei^^^ (8) 

q,a 

where Ska = — 2t(cos kx + cos ky) — ji — ah^ dj.{y) {k, q) = 
cos[A:^(^) - qx{y)/'A and 



k 

^l^'^ = ^T.d.iy){k,q){c-k+,lCkt). (10) 
k 

The resulting Hamiltonian can be diagonalized by the 
transformation 

Cka = (uknalna - Crvlnalla) , (H) 

n 

where jna are the quasiparticle operators while Ukna and 
Vkna fulfil the BdG equations 



(9) 



Ef J^ki 
\ A* 



Hkpa A^p 

X * rr* 

kp kpa 



^pncT 
^pna 



En. { )(.12) 

^km 



Here, 



J^kpa — (Jkp^ka ^k,p—QsDW r) ^'SjP+Qsdw ^ ' 



is the normal state Hamiltonian and 



(13) 



A _ J T.g^q^k,-p+q for S-WaVC, 

" I E« A" da{k, q)dk,-p+g for d-wave. ^'^> 
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Superconducting order parameters are determined self- 
consistently from Eqs. (|9ll0p together with 

n 

- ^U-k+q,niVl^^f{-Eni), (15) 
n 

where f{E) = [1 + exp(/3£^)]~^ is the Fermi-Dirac distri- 
bution function. In particular, we focus on the relation 
between the BCS A^^^^^ and FFLO A^qf^^^o) compo- 
nents of the superconducting OP. 

An unrestricted search for solutions of the BdG equa- 
tions in the momentum space has no advantage over the 
standard analysis in the real space. However, the ad- 
vantage becomes evident when the system is invariant 
under translations along one particular axis. Such a 
case will be considered in the present work, where we 
assume that the direction of FFLO modulation q is par- 
allel to incommensuration wave-vector of the magnetic 
order S = (7r,7r) — QsDW- Under this assumption Eq. 
([12]) represents eigenproblem of a block matrix. For the 
L X L system with 6\\d where d is a diagonal unit vector 
d = (1, 1), the blocks are not larger than 2Lx2L and each 
block account for particles and holes which momenta are 
connected by the following relation: p' = p + 27r dn/L^ 
where n = 0, — 1. In the second case, when 6\\x 
with X = (1,0), the biggest blocks consist of 4L x 4L 
elements. Within each block we consider particles and 
holes with momenta connected via p' = p + 27r xn/L^ as 
well as = p + (0, tt) -h 27r xn/L, where n = 0, L — 1. 

Certainly, we cannot exclude a possibility that an ad- 
ditional modulation in the perpendicular direction may 
lead to a further lowering of the free energy and stabiliza- 
tion of the superconducting phase beyond the boundaries 
obtained in the present studies. 

III. RESULTS AND DISCUSSION 

In the numerical calculations we take ja = — 0.2t, what 
gives the occupation number slightly below one electron 
per lattice site. The magnitudes of the pairing potentials 
are Vs = —2.0t for s-wave and Vd = —1.2t for d-wave 
pairing. The BdG equations have been solved for clusters 
up to 256x256 sites at temperature ksT = 10~H. 

The physically relevant solutions of the BdG equations 
are determined according to the following procedure: we 
start from a small system, e.g., 64 x 64 and for each mag- 
netic field h and SDW amplitude Mq we iteratively solve 
the BdG equations starting from LO states with all pos- 
sible momenta q. Depending on QsDW these momenta 
are either along (1,1) or (1,0) directions. Then, the so- 
lution with the lowest free energy together with the LO 
states with neighboring momenta are taken as the initial 
states in iterative solution of the BdG equations on much 
larger clusters, e.g., 256 x 256. Although our approach 
to the inter-site pairing remains valid for an arbitrary 



singlet state, we solve the BdG equations starting from 
an initial d-wave state with Af = — Aj^. 

Solving the BdG equations in the absence of SDW, 
we have found for the assumed dispersion relation and 
the Fermi energy that the FFLO state with vectors q 
along (1,0) direction have lower free energy than the 
states with q along (1,1). Therefore, we will show results 
only for SDW with the incommensuration wave-vector 
^11(1,0). However, we have found the same qualitative 
results also for 5||(1,1) and gr||(l,l), when the free en- 
ergy is only slightly higher. Note that we assume that 
the onset of incommensurate SDW does not change the 
direction of modulation of A^. It is the only restriction 
imposed on the solutions of the BdG equations in our ap- 
proach. Further on, when discussing the d-wave super- 
conductivity we present site-dependent superconducting 
order parameter defined for site i as an average of Aij 
over four bonds connecting site i to its neighbors: 

= ^ {^i,i-\-x + ^i,i-x — ^i,i-\-y ~ ^i,i-y) • (16) 
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FIG. 1. (Color online) Real-space profile of the supercon- 
ducting OP for QsDW = (yI , l)7r calculated on a 256 x 256 
cluster. Panels a and b show s-wave superconductivity with 
Mo = 0.2t for /i = and h = 0.25t, respectively. Panels c and 
d show the d-wave case with Mo = O.lt and h — (panel c) 
and h = 0.25t (panel d). 

Fig. [1] shows the real-space profiles of the supercon- 
ducting OP for s-wave (upper panels) and d-wave (lower 
panels) superconductivity in the absence (left column) as 
well as in the presence (right column) of the magnetic 
field. One can see that for = 0, A^ is spatially mod- 
ulated with periodicity given by the incommensuration 
6. In other words, considering translationally invariant 



BCS state in the presence of incommensurate SDW is an 
approximation and leads to a state with the free energy 
higher than that of the inhomogeneous states shown in 
Figs. [TJi and [it. From Eqs. (p fTQ|) one immediately 
finds that incommensurate SDW itself induces Cooper 
pairs with momenta along 6. The magnitude of SDW is 
expected to determine whether these components domi- 
nate or, as in the case shown in Fig. [TJ they represent 
corrections to the BCS superconductivity. 

Similarly to the case of non-magnetic systems, exter- 
nal field favors pairing with a non-zero momentum also 
in the presence of SDW. As demonstrated in the right 
panels in Fig. [H for sufficiently strong field, the BCS 
component Aq^o vanishes and the superconducting OP 
changes sign in the real-space. However, the spatial pro- 
file of is very different from a standard cosine depen- 
dence. Contrary to the LO phase, several components 
with different momenta of Cooper pairs give significant 
contribution to the superconducting OP. 




FIG. 2. (Color online) Phase diagrams for d-wave supercon- 
ductivity with QsDW = ( yI 5 1)^ obtained on a 64 x 64 clusters. 
Panel a shows the free energy F relative to the normal-state 
energy Fn under the same conditions SF — Fn — F. Panels 
b and c show the maximal max\Aq\ and the BCS |Ag| com- 
ponents of the superconducting OP, respectively, while panel 
d shows sum of all amplitudes with nonzero momentum of 
Cooper pairs J2q^o |A^|. 

Results shown in Fig. [1] indicate that incommensurate 
SDW affects the total momentum of Cooper pairs both 
in the presence and in the absence of external magnetic 
field. Then, the key question is whether these results are 
generic for incommensurate SDW or, in contrary, they 
require a fine-tuning of Mq and h. In order to answer 
this question we have calculated phase diagrams in the 
Mo-h plane. The diagrams shown in Figs. [2] and [3] are 
the main results of our study. In order to determine 
the properties of the superconducting state (especially 
the role of Cooper pairs with nonzero momentum) we 
have calculated the following quantities: SF = Fn — 
where F and Fn denote the free energy and the free 
energy of a nonsuperconducting (normal) state, respec- 
tively; max|Aq| - maximal single components of OP [see 
Eqs. (|9TTQl) ]: |Ao| = |Aq^o| - the BCS component of OP, 
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FIG. 3. (Color online) The same as in Fig. [2] but for s-wave 
symmetry (Ag Ag). 



and Xlq^o \ ^q \ ~ amplitudes with nonzero mo- 

mentum of Cooper pairs. Increasing the magnitude of 
incommensurate SDW causes reduction of |Ao| accom- 
panied by an increase of Xlg^o \ ^q\' "^^^^ holds true for 
both symmetries. It is a clear indications that incommen- 
surate SDW acts in detriment of BCS superconductivity 
and favors pairing with nonzero momentum of Cooper 
pairs. Since also max|Aq| decreases when Mq increases, 
this superconducting state is very different from the stan- 
dard LO superconductivity, when only two components 
with opposite momenta are relevant. 

A clear distinction between (i-wave (Fig. [2j) and s- 
wave (Fig. [3] ) superconductivity shows up in the phase 
diagrams when superconductivity is simultaneously af- 
fected by SDW and external magnetic field h. In the 
case of d-wave superconductivity there is no unusual in- 
terplay between these two mechanism and the role SDW 
is rather negligible as long as Mq is smaller than the 
critical field. This stands in strong contrast to the re- 
sults for 5-wave superconductivity. When considered 
separately, the external field and SDW are strong pair- 
breaking mechanisms. However, they are not so destruc- 
tive upon superconductivity when they emerge together. 
As we consider a real-space pairing with rather strong 
pairing potentials, our approach should be applicable 
to extremely type-II superconductors. Hence, it should 
be possible to explain the obtained results investigating 
the spatial variation of a local effective field defined by 
hef[{Ri) = h-\-Mo cos((5sDW^i) on a very short length- 
scale of the order of the coherence length. Let us define 
a fraction of sites 



'f' = ^T.(^{hc-\hMR^)\), (17) 
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and bonds 

= E ^ ('^^ - l^eff(i?,)|) e {he - |/leff(J?,)|) , 

{i,i> 

(18) 

where the effective field /leff is smaller than the critical 
field he determined in the absence of SDW. Here, 0{. . .) is 
the Heaviside step function. These quantities are shown 
in Fig. m One can see that at least in the case of s-wave 
superconductivity provides very simple explanation 
of the general structure of the phase diagram shown in 
Fig. [3l In the presence of incommensurate SDW, ex- 
ternal magnetic field h proportional to Mq increases the 
number of lattice sites where the effective magnetic field 
is smaller than he and in this way it effectively screens 
superconductivity against SDW. Of course, at the same 
time it increases the effective field /leff at some other sites, 
but in the case of the FFLO superconductivity the order 
parameter can be adjusted in such a way that the influ- 
ence of these sites is minimized. 
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FIG. 4. (Color online) Fraction of sites (a) and bonds (b), 
where the magnitude of the effective magnetic field is lower 
than he = 0.25t (a) and he = 0.29t (b). Results are obtained 
on a 128 x 128 lattice for Qsdw = (if , l)7r. 

In summary, we have solved the BdG equations for su- 



perconductivity coexisting with (assumed) incommensu- 
rate SDW. We have considered a case when spatial mod- 
ulations of the superconducting and antiferromagnetic 
orders takes place in the same direction. The transla- 
tional invariance in the perpendicular directions allowed 
us to study systems as large as 256 x 256. Our numerical 
data provide a clear evidence that incommensurate SDW 
itself diminishes the role of the BCS pairing and simul- 
taneously favors formation of Cooper pairs with nonzero 
total-momentum. It has recently been demonstrated that 
also a complementary dependence, where tendency to- 
ward formation of incommensurate SDW is enhanced by 
the presence of the FFLO order, is possible. These ob- 
servations hold true for s-wave and d-wave supercon- 
ductivity and supports the hypothesis that FFLO-type 
of superconductivity exist in the HFLT phase of heavy- 
fermion superconductor CeCoIns. 

For s-wave superconductivity we have found rather 
surprising result concerning the case when superconduc- 
tivity is affected simultaneously by external magnetic 
field and the incommensurate SDW. We have demon- 
strated that these two mechanisms are less destructive 
upon s-wave superconductivity when they emerge to- 
gether. It means that under external magnetic field su- 
perconductivity may coexists with stronger SDW than in 
the absence of magnetic field. 
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